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We report a new type of phononic crystals with topologically non-trivial bandgaps for both longit¬ 
udinal and transverse polarizations, resulting in protected one-way elastic edge waves. In our design, 
gyroscopic inertial effects are used to break the time-reversal symmetry and realize the phononic 
analogue of the electronic quantum Hall effect. We investigate the response of both hexagonal and 
square gyroscopic lattices and observe bulk Chern number of 1 and 2, indicating that these struc¬ 
tures support single and multi-mode edge elastic waves immune to back-scattering. These robust 
one-way phononic waveguides could potentially lead to the design of a novel class of surface wave 
devices that are widely used in electronics, telecommunication and acoustic imaging. 


Topological states in electronic materials, including 
the quantum Hall effect pQ and topological insulat¬ 
ors HE], have inspired a number of recent develop¬ 
ments in photonics HE], phononics |6HIQ] and mechan¬ 
ical metamaterials mm- In particular, in analogy to 
the quantum anomalous Hall effect m , one-way wave¬ 
guides in two-dimensional photonic systems have been 
realized by breaking time-reversal symmetry 1T6UT8] . 
Moreover, very recently unidirectional edge channels 
have also been proposed for scalar acoustic waves by 
rotating fluids [19, 20 and for elastic waves by Coriolis 
force in a non-inertial reference frame m • However, 
the latter system is very difficult to implement in solid 
state devices. 

Here, we introduce gyroscopic phononic crystals, 
where each lattice site is coupled with a spinning gyro¬ 
scope that breaks time-reversal symmetry in a well- 
controlled manner. In both hexagonal and square lat¬ 
tices, gyroscopic coupling opens bandgaps character¬ 
ized by Chern numbers of 1 and 2. As a result, at 
the edge of these lattices both single-mode and multi- 
mode one-way elastic waves are observed to propagate 
around arbitrary defects without backscattering. 

To start, we consider a hexagonal phononic crys¬ 
tal with equal masses (m 2 = mi) connected by lin¬ 
ear springs (green and black rods in Figs. 0- Of 
and Eh The resulting unit cell has four degrees of 
freedom specified by the displacements of mi and m 2 
(U = [uff 1 , u™ 1 , u ™ 2 , u™ 2 }). Consequently, there are 
a total of four bands in the band structure (Fig. Eh 
Note that this is the minimal number of bands required 
to open a complete bandgap, since the first two elastic 
dispersions are pinned at zero frequency. The phononic 
band structures are calculated by solving the dispersion 
equation |22] 

[K(n) - w 2 M] U = 0 (1) 

for wave vectors pi within the first Brillioun zone. Here, 
uj denotes the angular frequency of the propagating 
wave and M = diagjmi, mi, m 2 , m 2 } is the mass 
matrix. Moreover, K is the 4x4 stiffness matrix as a 


function of the Bloch wave vector /x. The band struc¬ 
ture of this simple lattice is shown in Fig. [l^ and has 
a quadratic degeneracy between the third and fourth 
bands at the center of the Brillouin zone and a com¬ 
plete gap between the second and third bands due to 
the lack of inversion symmetry. However, this gap is 
topologically trivial, since time-reversal symmetry is 
not broken and the Chern numbers of the bands are 
all zero. 

In order to break time-reversal symmetry, we intro¬ 
duce gyroscopic coupling [23H25] and attach each mass 
in the lattice to the tip of the rotational axis of a gyro¬ 
scope, as shown in Fig. El- Note that the other tip 
of the gyroscope is pinned to the ground to prevent 
any translational motion, while allowing for free rota¬ 
tions. Because of the small-amplitude in-plane waves 
propagating in the phononic lattice, the magnitude of 
the tip displacement of the each gyroscope is given by 

Utip = hsinO « hO = hQe luJt for |0| <C 1, (2) 

where h and 6 denote the height and nutation angle 
of the gyroscope (Fig. and 0 is the amplitude of 
the harmonic change in 0. Interestingly, the coupling 
between the masses in the lattices and the gyroscopes 
induces an in-plane gyroscopic inertial force perpendic¬ 
ular to the direction of U tiv [24] [26J [27 : 

F g = ±iuj 2 aU Up , (3) 


where a is the spinner constant that characterizes the 
strength of the rotational coupling between two inde¬ 
pendent inertias in the 2D plane. As a result, the mass 
matrix in Eq. 0 becomes, 


M = M + 


0 

—ia 1 


0 

0 


ia.\ 0 0 \ 

0 0 0 

0 0 IOL2 

0 —ia 2 0 ) 


( 4 ) 


where aq and a 2 denote the spinner constants of the 
gyroscopes attached to m\ and m 2 , respectively. We 
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Figure 1: Ordinary and Gyroscopic Phononic Crys¬ 
tals: (a) Schematic of the hexagonal lattice. The purple 
and grey sphere represent concentrated masses m\ and m 2 , 
respectively. The green and black straight rods represent 
mass-less linear springs with stiffness k\ and k<i — fci/20, 
respectively. The dashed cell is the primitive cell of the lat¬ 
tice. (b) Schematic of a gyroscope with the top tip pinned 
to a mass in the lattice, (c) Unit cell for the ordinary 
(non-gyroscopic) phononic crystal, (d) Unit cell for the 
gyroscopic phononic crystal, (e) Band structure of the or¬ 
dinary (non-gyroscopic) phononic crystal (an = <22 = 0). 
The inset is the Brillouin zone, (f) Band structure of the 
gyroscopic phononic crystal (an — ol^ — 0.3mi) with the 
Chern numbers labeled on the bulk bands. The frequencies 
are normalized by ujq — \fk^Jm\. 


note that the imaginary nature of the gyroscopic iner¬ 
tial effect indicates directional phase shifts with respect 
to the tip displacements, which breaks time-reversal 
symmetry. 

We now consider a gyroscopic hexagonal lattice with 
<^2 = on = 0.3mi and show its band structure in Fig. 
[if. Comparing this to the original band structure re¬ 
ported in Fig. [J, we observe that the original quad¬ 
ratic degeneracy, between the third and fourth bands, 
is opened into a full band gap. Moreover, we also 
find that the original bandgap between the second and 
third bands first closes and then reopens as we gradu¬ 
ally increase the magnitude of oq and a^. In partic¬ 
ular, for oi\ = c ^2 = 0.07mi the gap is closed and a 
pair of Dirac cones at K points emerges (see Fig. SI 



Figure 2: Edge modes in Gyroscopic Phononic Crys¬ 
tal: (a) ID band structure showing bulk bands (black dots) 
and edge bands (colored lines). Red solid lines represent 
edge modes bound to the top boundary, while blue dashed 
lines represent edge modes bound to the bottom bound¬ 
ary. (b) Modal displacement fields of top edge states with 
negative group velocities, (c) Modal displacement fields of 
bottom edge states with positive group velocities 


for details). Importantly, for 0 L 2 = aq = 0.3mi (Fig. 
[if) both bandgaps are topologically-nontrivial as high¬ 
lighted by the non-zero Chern numbers labeled on the 
bands (the calculations conducted to compute these to¬ 
pological invariants are detailed in the Supplementary 
Materials E3). Therefore, in the frequency ranges of 
these nontrivial bandgaps, we expect gapless one-way 
edge states, whose number is dictated by the sum of 
Chern numbers below the bandgap. 

To verify the existence of such one-way edge states, 
we perform one-dimensional (ID) Bloch wave analyses 
on a supercell comprising 20 x 1 unit cells, assuming 
free boundary conditions for the top and bottom edges. 
In full agreement with the bulk Chern numbers, the 
band structure of the supercell shows one one-way 
edge mode on each edge in both bandgaps. For modes 
bound to the top edge (Fig. If)- the propagation 
can only assume negative group velocities (red solid 
lines with negative slope in Fig. If)- On the other 
hand, the modes bound to the bottom edge (Fig. [2f) 
possess positive group velocities (blue dashed lines 
with positive slope in Fig. [2fi). Since these edge modes 
are in the gap frequency range where no bulk modes 
may exist, they cannot scatter into the bulk of the 
phononic crystal. Furthermore, their uni-directional 
group velocities guarantee the absence of any back 
scattering and result in the topologically protected 
one-way propagation of vibration energy. 

To show the robustness of these edge states, we con¬ 
duct transient analysis on a finite sample comprising 
20 x 20 unit cells with a line defect on the right bound¬ 
ary created by removing twelve masses and the springs 
connected to them (Fig. If)- A harmonic force excit¬ 
ation, Foe - ^, is prescribed at a mass site on the bot¬ 
tom boundary (red arrow in Fig. [3fi) with frequency 
within the bulk bandgap between the second and third 
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Figure 3: Transient Response of a gyroscopic phononic 
crystal consisting of 20 x 20 unit cells with a line defect on 
the right boundary: Snapshots of the displacement field at 
(a) t = 2To, (b) t = 12To, (c) t = 22To and (d) t = 32To, 
where To = is the characteristic time scale of the 

system. Starting from t — 0, a time-harmonic excitation 
force F (t) = [F x (t) , F y (t)\ = [1, l]Foe~ ZUJt is prescribed at 
the site indicated by the red arrow. 


bands (uj/uj o = 1.3). In Fig. [3] we plot snapshots of 
the velocity field at different time instances, t/T 0 = 2, 
12, 22 and 32, where To = yfmi/ki is the character¬ 
istic time scale of the system. Remarkably, because of 
their topological protection, the edge modes circum¬ 
vent both the sharp corner and the line defect without 
any reflection. 

Next, we investigate the effect of the lattice geometry 
and start with an ordinary square phononic crystal 
with masses mi connected by springs with elastic con¬ 
stant k\. To make the lattice statically stable, we add 
an additional mass m 2 = mi at the center of each unit 
cell and connect it to its four adjacent mi masses by 
springs with elastic constant k^ = 2£q (See Figs. & 
andEb). The band structure for this lattice (shown in 
Fig.lml) contains a pair of three-fold linear degeneracy 
among the first, second and third bands at the X points 
of the Brillouin zone. Note that this type of degener¬ 
acy, consisting of a locally flat band and a Dirac point, 
is known as the “accidental Dirac point” [28]. Interest¬ 
ingly, while previously this degeneracy has been found 
to be very sensitive to the system parameters [28], in 
our lattice it robustly appears at the X points when 
mi = m 2 . Upon the introduction of gyroscopic inertial 
effects (c ^2 = a 1 = 0.3mi), these three-fold degener¬ 
ate points are lifted and a gap is created between the 
second and third bands. The Chern numbers of the 
two bulk bands below the gap is two, predicting the 
existence of two topological edge states. In Figs. [4]f 
and [4jy, we plot the band structure of the correspond¬ 
ing 20x1 supercell, highlighting the two one-way edge 
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Figure 4: Square Lattice results: (a) Schematic of the 
square lattice. The purple and grey sphere represent con¬ 
centrated masses mi and m 2 , respectively. The black and 
green straight rods represent mass-less linear springs with 
stiffness k\ and fe = 2 k±, respectively, (b) Unit cell 
for the ordinary (non-gyroscopic) phononic crystal, (c) 
Unit cell for the gyroscopic phononic crystal, (d) Band 
structure of the ordinary (non-gyroscopic) phononic crys¬ 
tal (an = 02 = 0). The insets is the Brillouin zone, (e) 
Band structure of the gyroscopic phononic crystal (cq = 
a 2 = 0.3mi) with the Chern numbers labeled on the bulk 
bands. Frequencies are normalized by cjo — y/ki/rrii. (f) 
ID band structure showing bulk bands (black dots) and to¬ 
pological edge bands (red solid lines for edge modes bound 
to the top boundary. Note that the edge bands that are 
bound to the bottom boundary are not shown here, (g) 
Modal displacement fields of the edge states shown in (f). 


modes and their modal displacement fields. 

To summarize, we demonstrated that gyroscopic 
phononic crystals can support topologically non-trivial 
gaps, within which the edge states are unidirectional 
and immune to back-scattering. Moreover, transient 
analysis confirmed that the propagation of such topo¬ 
logical edge waves is robust against large defect and 
sharp corners. We note that phononic crystals [29 , [30] 
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and acoustic metamaterials [3X1(35] that enable manip¬ 
ulation and control of elastic waves have received signi¬ 
ficant interest in recent years [221(36] . not only because 
of their rich physics, but also for their broad range 
of applications [371(49] . Interestingly, the edge wave 
modes in phononic crystals are important in many 
scenarios [501152] . including vibration control [53] and 
acoustic imaging [45]. However, most of reported stud¬ 
ies have focused on topologically trivial surface waves 
that can be easily scattered or localized by defects [50J. 
Therefore, the work reported here could open new av¬ 
enues for the design of phononic devices with special 
properties and functionalities on edges, surfaces and 
interfaces. 
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